Nonexistence of intrinsic spin currents 



o 
o 

wo: 
< 



i 

C 

o 
o 



> 
m 
oo 

o 

O 

-i— > 



i 

C 

o 
o 



Alexander Khaetskii 
Physics Department and Center for NanoScience, LMU Miinchen, D-80333, Germany 

(Dated: February 2, 2008) 

We have described the electron spin dynamics in the presence of the spin-orbit interaction and 
disorder using the spin-density matrix method. We showed that in the Born approximation in the 
scattering amplitude the spin current is zero for an arbitrary ratio of the spin-orbit splitting and 
the scattering rate. Various types of the disorder potential are studied. We argue that the bulk spin 
current has always an extrinsic nature and depends explicitely on scattering by impurities since it 
appears only beyond the Born approximation in the scattering amplitude. 

PACS numbers: 72.25.-b, 73.23.-b, 73.50.Bk 



Spin-orbit coupling brings about a number of interest- 
ing effects, one of which is generation of a spin flux in the 
plane perpendicular to the charge current direction. This 
phenomenon occurs in the paramagnetic system and is 
very well known for quite a long time, see Ref. 1], where 
the Yafet-Elliott spin-orbit mechanism was considered. 
It is a consequence of the fact that in the presense of 
spin-orbit coupling the scattering by impurities has an 
asymmetric character (the Mott effect) [2|. Spins with 
up-orientation are scattered preferably to the right and 
with down-orientation - to the left. This phenomenon 
exists only beyond the Born approximation in the scat- 
tering amplitude and leads to an accumulation of the 
spin density near the sample surface Q- Mutual trans- 
formation of the current and spin fluxes leads also to the 
renormalization of the electrical conductivity of the sys- 
tem, see Ref. 0] , where the case of 3D holes described by 
the Luttinger Hamiltonian was considered. It has been 
recently claimed 0, that an analogous phenomenon 
can exist even without scattering by impurities, i.e. in 
the ballistic regime, the corresponding contribution being 
called "intrinsic" or dissipationless. Later several papers 
appeared where the effect of scattering by impurities was 
taken into account 0, S @ with a range of totally dif- 
ferent results. This was done by the Kubo formula. We 
solve this problem using the well known method of a spin- 
density matrix |3j. We argue that intrinsic spin current 
in the bulk cannot exist. Any spin current must be due 
to asymmetric scattering by impurities (the Mott effect). 
Since this phenomenon occurs only beyond the Born ap- 
proximation in the scattering amplitude, this contribu- 
tion is explicitely dependent on the impurity scattering 
and leads to the well known extrinsic contribution. 
In particular, there cannot exist the universal value of the 



spin current since it does not depend on the scattering 
properties. In the Born approximation (when the scat- 
tering amplitude has additional symmetry propeties, see 
below) the Mott effect is absent, and the spin current is 
zero for an arbitrary value of At, where A is the spin 
splitting of the electron spectrum and r is the transport 
scattering time. We have shown this by exact calcula- 
tions for the case of the Rashba Hamiltonian. [ljj Thus 
the correct terminology, i.e. intrinsic or extrinsic contri- 
bution, should be used in accordance with the strength of 
the scattering (Born or beyond) rather then the presence 
or absence of the impurity scattering itself. Therefore, 
certain care should be taken when trying to check by 
"exact" numerical diagonalization the robustness of the 
" intrinsic" value with respect to the disorder since strong 
scattering inevitably generates an extrinsic contribution 
to the spin current. 

The Hamiltonian of the problem is 

2 2 

W(P) = + e M (p) = ^- + Map, (1) 

2m 2 2m 

where fi(p) = [n • p], n is the unit vector normal to the 
2D plane (z-axis), cm(p) are the eigenvalues, M = ±1/2 
are the helicity values. The eigenfunctions are 



lU=±l/2 |U=±l/2 

where D^p(Sl) is the rotation matrix 0], <t> the angle of 
p, and the eigenfunction of the a z operator. 

Spin current, kinetic equation. We will calculate the 
q yz component of the spin current. This quantity is zero 
in the thermodynamic limit and defined as 



qyz =TrJ AjL f ( p )l(S z V y + V y S z ) =TrJ AJLbLf(p)§, cx (/+_ + /_+). (2) 



r 



Here /(p) is the spin density matrix 12], V y the y- component of the velocity operator and S z = (l/2)a z 
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the spin operator. The last expression in Eq.JSJ is given 
in the helicity basis. The general expression for the quan- 
tum kinetic equation in the case of spin-orbit interaction, 
when the Hamiltonian and the Wigner distribution func- 
tion are matrices over the spin indexes, was derived in 
Ref. 3] . When there is a magnetic field or some inhomo- 
geneity in the problem, the field term and the gradient 
term must be symmetrized since the velocity operator is 
also a matrix. In our case when we deal only with the 
electric field which is constant in space this equation is 
simple and reads 



the expression for the collision term is given below. Now 
we write Eq.|0 in the helicity basis where the Hamilto- 
nian is diagonal. While doing that, we should take into 
account the fact that eigenfunctions XMp depend on the 
direction of the momentum p, thus the matrix elements 
of the derivative df /dp in this basis do not coincide with 
the quantities <9/mm'/^P 



9/mm' i r ?i . % * d\M'p 

— Tl a > j\mm>, a-MW =inx Mp - 




MM' 



dp 



dp 



+ eE !p~ + i [H to> ft = st{I{p)} (3) 

The last term on the left hand side is a commutator and 



We see that there appears the commutator of the vector 
matrix a with /. Thus for Eq.@ in the linear response 
regime (E || x) we obtain 



, ; v ( , s ^ MM , _ i s J^ eE (fW M ,(p) - /<°k (p)) + Ue M {p) - e M ,(p))f MM ,(p) = St{f(p)) MM , 



dp 



P 



(4) 



Here f^[ M {p) is the equilibrium Fermi function corre- derived in many papers, for the refs. see and in 

sponding to the helicity value M. The collision term was the helicity basis has the form 



st(Rp)) 



MM' 



d 2 p 



{2nh) 



M\M[ 



-6(e Ml (p) ~ e Mi ( Pl ))[K^,(u; ppi ) ■ f Ml M'(p) + JmmAp) • *$$(<<w)]}> (5) 
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where the kernel in the Born approximation in the scat- 
tering amplitude is: 



The additional symmetry property indicated here exists 
only in the Born approximation. [2| 



,(1/2)* 



, ., 7,- 2 (1/2)* Smooth scattering potential. First consider the math- 

K P1 ) = j.\u{p- "Pi)l m ND MMi (Wp P1 )Z> M , M , (w PPl )iematically simple case of a smooth scattering potential 

(6) when the intcrband transitions are suppressed which is 
realized at maR/h 3> 1, R being the radius of impurity. 
Then from Eas. (|4l5|l we obtain 



Here N is the 2D impurity density, U(p — pi) 
is the Fourier component of the impurity potential. 



(1/2 

D ( M ' M ' (w PPl ) depends only on the scattering angle 9 



Diagonal components 



(1/2) 
/2.1/2 



D 



(1/2) 
^1/2,-1/2 



cos(0/2), and D { ^l l/2 = £» < L 1 / /2.i/2 = -» sin(6>/2). The 



1/2,-1/2 -1/2,1/2 

Born scattering amplitude is given by 



< p P1 «^(«^pxMp-pi)i K^-(FM P pi r- 



,(1/2) , 



(7) 



eE 



eE 



dfl_ 
dp 



df. 



(0) 



dp 



^/++ + ^(/+- -/-+), (8) 
v+ v + 

£(/+-/-+), (9) 



i eE 
Yp~ 



{Si 



/°) + ^(e+-e_)/ + _ 



c(^/++ 



(10) 
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vJ 



(ii) 



where d = a,b = —c, c = —ia, a = 
-±Jd9/(2n)W(9)sm 2 9, quantity -2ap/V+ is 
equal to the inverse transport scattering time r, 
W{6) = N ■ \U(p - Pl )\ 2 /2h 3 . #(p),/£(p) are the 
equilibrium Fermi functions which correspond to the 
helicity ±, V±(p) = p/m ± a/2 are the velocity values 
for a given p for ± bands. The expressions for the 
coefficients a,b,c,d are exact but should be used here 
for 9 <C 1 since we consider small- angle scattering. In 
Eqs.JHEU th e quantities f++(p), f+-(p), f-+(p), f—{P) 
depend only on the modulus of p. In deriving these 



Eqs. we used the following angular dependences of the 
components of the matrix /(p): /++(p), /--(p) oc cos0 

and /_| (p),/ ^(p) sin</>. Besides, we used the 

symmetry properties of the matrix K(9): 



A +4 



K~ 



(12) 



which can be easily proved from Eas. (|6l7|l . The quan- 
tities entering Eas. (|8lll() have the following relations to 
the average spin components: 



J 



< S z >cx (/+_ + /_+), < S ■ p >tx (/+_ - /_+), < S ■ n >oc (/++ - /_ 

I 



(13) 



The last quantity exists even in the thermodynamic limit 
(for a given momentum p) . From the above equations for 



the quantity of interest we find 



J 



dp dp 



"All -f-) = -ke+ -e_)(/+_ +/_+). 



(14) 



r 



This equation is exact for an arbitrary value of At, 
A = (e + — e_) = ap. Hence, q yz = 0, i.e. the spin 
current is zero. Note that Ea.l|14fl has a clear phys- 
ical meaning. The second term on the left hand side 
was taken into account before and describes the ap- 
pearence in the electric field of the z-component of the 
spin due to the angular dependence of the wave functions. 
Exactly this term gives the contribution e/8ir after inte- 
gration in Eq.J2J). However, the first term in Eq. i|14|) 
describes the change in the distribution functions due to 
the accelaration along the electric field and cancels ex- 
actly the contribution of the second term after integra- 
tion in Eq.©. Q Note that in Eas.l|8l9fl the scattering 

admixes only the component (/_| — / |_). This is the 

direct consequence of the Born approximation and the 
symmetry properties, Eci. (|12f) . Beyond the Born approx- 
imation the quantity + ^++ ) is n °t zero. Exactly 
this quantity is responsible for the generation of the spin 
flux due to the scattering when the particle flux flows 
in the sample, see Ref.[3j. When this quantity is not 
zero, in Eqs. (US© there appears the term proportional to 
(/-I + / |_) which means the appearence of q yz due to 



the Mott effect when the current flows in the x-direction. 

5-scattering potential. Here we consider the case of a 
short range scattering potential when W{9) = Wo (con- 
stant). Then interband transitions are allowed and from 
Eas. 1)415(1 we obtain a system of coupled equations similar 
to Eas. (|8lll|) where now the components of spin-density 
matrix for the values of p± = p ± ma appear (see Fig.l). 
To simplify the presentation we will consider only the 
limiting cases of large and small At. When Ti/t <C A 
we find that (/+_ - /_+)/(/+_ + /-+) * (Ti/t)/ A <C 1. 

Neglecting everywhere the (/_| — / |_) components, for 

the z-component of the spin we obtain 



— (f + (p) ~ f-(p)) + Ue + (p) 
p n 



e _(p))(/+_(p)+/_ + (p)) 



2cXr, 



2cY p , 



where 

x P = ^-^M, y^^u+iP-y^f--^) 

(15) 

Note that in Ea. H5fl the velocities V±(p) enter at the 
momentum p since for horizontal transition the velocity 
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FIG. 1: Schematics of the ± energy bands. Momenta p,p±,Pp 
are shown, see the text. 



is conserved. For the /++,./- 
equations 



components we have the 



j? < \ 2a P t t \ , 2a iP+ 
E+(p) = —f++(p) + — rr— 
v + v+ 



(/__(?+) + 2/++(p)), 



(/ ++ (p_)+2/__(p)) J 



where we introduced the notations: E+(p) — 
eEdfl(p)/dp, £L(p) = eEdf°(p)/dp and again V±(p) 
enter at the momentum p. ai = ^j 1 J d#/(27r)(l — 
cos 9) cos = — Wo/4. Writing these Eqs. for the 
momenta values p±, using cii = a and the relations 
V+(p) = K_(p + ), F+(p_) = /°(p) = /°(p+), 

/+(?-) = /-(p)> we obtain 

P+f++(p) =Pf — (P+), P-f—ip) =Pf++(P-)- (16) 
Then the solution is 

,PV+ 1 



f++(p) = e +(pY 



2a {p + p+) 



2 ■ 



(17) 



and the expression for / (p) is obtained from Eq. i|17|) 

by replacing + by — . Using Eas. (|17ll6|) we can calculate 
quantities X p ,Y p and with the use of c = — ia finally 
obtain 



p \2p + ma 2p — ma 



--(£+-€_)(/+_+/_+) 



(18) 



Here the ratio between ma and p is arbitrary, the only 
restriction is that the position of the Fermi level should 
correspond to the values of p where the above men- 
tioned equalities between the velocities are still valid, 
see Fig.l. Integrating in Eq.© between the points 
Pp = + \/pp + (ma) 2 /A, see Fig.l, we again ob- 

tain q yz = 0. In the opposite case A = we immediately 

see from Eas. (|4l5|l that (/+_ + / h) = and the spin 

current is zero. 

In conclusion, using the spin-density matrix method 
for the case of the Rashba Hamiltonian we showed that 
within the Born approximation in the scattering ampli- 
tude the intrinsic spin current is zero for an arbitrary 
ratio of spin splitting and the impurity scattering rate. 
We argue that the spin current appears only beyond the 
Born approximation, depends explicitcly on the scatter- 
ing and corresponds to the well known extrinsic spin cur- 
rents [10. 

After this work had been completed I became aware of 
recent work [l^J . It is not clear to me to what extent the 
authors generalize their conclusion about the absence of 
the spin currents in the bulk. Again, my opinion is that 
extrinsic dc spin currents can flow even in the bulk P, Q . 

I am grateful to M.I. D'yakonov, L. Glazman and E.I. 
Rashba for fruitful discussions and also to the partici- 
pants of the PASPSIII conference for the interest to my 



work. 
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